Elliptic flow (^2) in pp collisions at energies available at the CERN Large Hadron 

Collider: A hydro dynamical approach 



S. K. Prasad, Victor Roy, S. Chattopadhyay, and A. K. Chaudhuri 
Variable Energy Cyclotron Centre, 1/AF Bidhan Nagar, Kolkata 700 064, India 
(Dated: September 28, 2010) 

At Large Hadron Collider energy, the expected large multiplicities suggest the presence of collec- 
tive behavior even in pp collisions. A hydrodynamical approach has been applied to estimate the 
expected elliptic flow measured by the azimuthal asymmetry parameter V2 in pp collisions at y^s = 
14 TeV. V2 of 7r~ is found to be strongly dependent on the parton density profile inside a proton 
[e.g., surface diffuseness parameter (5)]. For ^ = 0.105 fm, V2 is found to be positive and larger 
compared to that at ^ = 0.25 fm. The centrality dependence of V2 has also been studied. 
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I. INTRODUCTION 

The elliptic flow is found to be one of the most inter- 
esting observables at the relativistic heavy-ion collider 
(RHIC). Large elliptic flow in noncentral Au-|-Au colli- 
sions confirms fluidlike behavior of the produced mat- 
ter [ij . Elliptic flow measures the momentum anisotropy 
of the produced particles. In noncentral collisions be- 
tween two identical nuclei, the reaction zone is spatially 
asymmetric. The rescattering process among the pro- 
duced particles (locally isotropic in momentum space) 
transfers this spatial asymmetry into the momentum 
space, and the momentum distribution of the produced 
particles becomes anisotropic. Elliptic flow is an early 
time phenomenon. It is a sensitive probe to the: (i) de- 
gree of thermalization, (ii) equation of state (EOS), and 
(iii) transport coefficients [2|t3- 

In pp collisions at RHIC energy, the average multiplic- 
ity (rich) is significantly low for any collective effects to 
be seen, however, detailed studies of various observables 
e.g., Hanbury-Brown-Twiss parameters Q and charged 
particle spectra gave hints that even in pp collisions, col- 
lective models (e.g., hydrodynamics) might give a sat- 
isfactory explanation for a large part of the data. Ap- 
plicability of hydrodynamics in a small system such as 
pp is uncertain. Hydrodynamics requires local thermal 
equilibration, which can be achieved only if the mean- 
free path of the constituents is small compared to the 
size of the system A << R. In pp collisions, size of the 
system is not large A ~ R ~ 1 fm. However, in Ref. 
it is argued that, if the medium is isotropized within a 
time scale n, hydrodynamics may be applicable beyond 
Ti. See Ref. [1| and references therein, where the appli- 
cability of hydrodynamics in central pp collisions after 
~ 0.2 fm is justified. The scenario of collective expan- 
sion of matter created in pp collisions at Larce Hadron 
Colhder (LHC) energy is discussed in Ref. Q- The ex- 
pectation that hydrodynamics can be applicable in high- 
energy proton-proton collision is not new, for example, 
see Ref. Recently, the possibility to observe a col- 

lective expansion signal - in the form of an azimuthal 
anisotropy of particle production with respect to the re- 
action plane - caused by multipartonic interactions in 
proton-proton collisions at the LHC is studied [ll|. In 



Ref. ^ll'l, different values of integrated V2 are predicted, 
which range from -3% to 10% depending upon the profile 
of matter distribution in transverse space of the collid- 
ing protons, and it is argued that the study of hadron 
anisotropy with respect to the reaction plane in pp colli- 
sions at LHC energies can provide important information 
on the proton shape and structure. In Ref. [l3|, large el- 
liptic flow ~10-20% is predicted in high-multiplicity 
pp collisions at the LHC if hot-spot-like structures are 
produced in the initial collisions. 

In pp collisions at LHC energies, average multiplicity 
will increase, and there could be events with multiplic- 
ity comparable to the multiplicity in peripheral Au-|-Au 
collisions at the RHIC. A picture based on two separate 
transverse distance scales in pp collisions at higher ener- 
gies gives an impact-parameter dependence of inelastic 
(INEL) collisions. In the case of central pp collisions, the 
distribution of hard partons {x > 10^^) in the two collid- 
ing nucleons will overlap, while in large impact-parameter 
collisions, partons with x << 10^^ will overlap with sig- 
nificant probability |13j]. A trigger on hard dijet pro- 
duction can quantitatively distinguish between a central 
and a peripheral collision. In this scenario, these high- 
energy pp collisions can pictorially be represented as AA 
collisions, and the impact-parameter dependence of ini- 
tial eccentricity Sx can be estimated. Significantly large 
energy density, which generates large multiplicity in pp 
collisions at the LHC will lead to rescatterings, thereby 
creating a system where hydrodynamics can be applied 
and the parameters such as V2 can be estimated. Here 
we report the results obtained by applying the hydrody- 
namic evolution to the system formed in pp collisions at 
the LHC energy (^/s = 14 TeV). 

The paper is organized in the following way. Section II 
describes the details of the code (AZYHYDRO) used for 
the hydrodynamic evolution of the system, and we ex- 
plain how we incorporated the impact-parameter depen- 
dence of the pp INEL cross section. In Sec. Ill, the results 
are presented and are discussed. We have explained how 
the initial energy density is chosen for the hydrodynamic 
evolution of the system at -^/s = 14 TeV. The pt spectra, 
mean pt, dN/dy, spatial eccentricity, and V2 with their 
dependence on impact parameter b or pt are presented. 
In Sec. IV, we summarize our results and conclude. 
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II. THE HYDRODYNAMIC MODEL 

We briefly describe the hydrodynamic model used 
to obtain space-time evolution of fluid formed in pp 
collisions. Details can be found in Refs. EH. The 
equation of motion of a relativistic ideal fluid follows 
from the local conservation laws of energy and momen- 
tum, and other conserved currents (e.g., baryon number), 



d^T^^ix) = 
d^f{x) = 0. 



(1) 
(2) 



Ideal fluid decompositions of energy-momentum tensor 
{T^^) and baryon four-current (j^) arc as follows: 

Tt'^ix) = [e{x)+p{x)]u^'{x)u''{x)- g^'^'pix), (3) 
f{x) = n{x)ui'{x), (4) 

where e{x) is the energy density, p{x) is the pressure, 
and n{x) is the conserved baryon number density at 
point x^^ = {t,x,y, z}. is the hydrodynamic four- 
velocity, = 7(l,W:r,Wy,W2) with 7 = ^(^_^/_^2„^2) ' 

The publicly available AZYHYDRO code solves" 2-fl- 
dimensional hydrodynamics. At high collision energies, 
relativistic kinematics and its influence on the particle 
production process implies longitudinal-boost invariance 
of the collision fireball near midrapidity jl^. As a 
result, the longitudinal velocity field scales as z;^ = |, 
and it is convenient to use a coordinate system spanned 
by longitudinal proper time r = t^l — v1 and the 
space-time rapidity rj = instead of t and z. 

Longitudinal-boost invariance is then equivalent to 
77 independence. By assuming the longitudinal-boost 
invariance, we reduce the number of energy-momentum 
conservation equations from four to three, viz. two 
transverse and one time. This will restrict us as we 
see the effect of rapidity dependence of transverse flow 
pattern. 

The initial thermalization stage lies outside the domain 
of applicability of hydrodynamical approach and must be 
replaced by the initial condition for the hydrodynamical 
evolution. We assume that the fluid is thermalized at 
the initial time r^. We have systematically performed 
the calculations for three different values of (0.2, 0.4, 
and 0.6 fm). At the initial time, transverse velocity of 
the fluid is zero Vx{x,y) = Vy{x,y) = 0. At an impact- 
parameter b, the initial energy density is assumed to be 
distributed as 



y, b) = eoNcouix, y, b). 



(5) 



where Ncoii{x, y, b) oc T{x + ^,y)T{x— ^,y) is a Glauber 
model calculation for the transverse profile of the (par- 
tonic) binary collisions at impact parameter b. The cen- 
tral energy density eg does not depend on the impact 
parameter of the collisions. 



For the Glauber model calculation of collision number 
distribution, we assume a Woods-Saxon profile for the 
(partonic) density distribution of the colliding protons. 



p{r) 



Po 



1 + e('-«)/« 



(6) 



with R=1.05 fm. The calculations are performed at two 
values of proton diffuseness parameter (i.e., ^ = 0.105 fm 
and ^ = 0.25 fm) (see Ref. and references therein) . 
The value of po is fixed in such a way that the total INEL 
cross section in pp collisions at a given energy is repro- 
duced in the Glauber model calculation. The values of 
Po, required to produce the INEL cross section are 0.45 
and 0.70 at -^/s = 0.2 and 14 TeV, respectively, for ^ = 
0.25. The total INEL cross sections in pp collisions at y/s 
— 0.2 and 14 TeV are taken to be 40 and 80 mb, respec- 
tively. The parton-proton thickness function is given by 
the optical path length: 



T{x,y) 



p{x, y, z)dz. 



(7) 



Hydrodynamical equations EqlU and [2] are closed only 
with an EOS p — p{e,p). We have used the EOS, which 
is composed of lattice EOS and hadron resonance gas 
EOS. Recently Cheng et.al. [l7| presented high statistics 
lattice QCD results for the bulk thermodynamic observ- 
ables (e.g., pressure, energy density, entropy, etc.). We 
have parametrized the entropy density as 



(8) 



The values of the parameters a, /3, 7 are chosen in such a 
way that the lattice simulation of s/T^ is best fitted fl8| . 
We have taken the crossover temperature = 196 MeV. 
From the parametric form of entropy density, pressure 
and energy density can be obtained by using the thermo- 
dynamic relations: 



p{T)^ Js{T')dT', 
e(T) = Ts-p. 



(9) 
(10) 



We complement the lattice EOS [17l| by a hadronic res- 
onance gas (HRG) EOS, which comprises all the reso- 
nances below the mass 2.5 GeV. The entropy density of 
complete EOS is obtained as. 



with X — ^-5^ 



s — 0.5[1 -t- tanh{x)]sHiiG + 0.5[1 — tanh{x)]s lattice 

(11) 

= 0.1 Tc. shrg is the entropy density 
for HRG, and slattice is the entropy density for lattice. 

As the system expands, its volume increases, and the 
density decreases, as a result, after some time, the mean- 
free path of the particles becomes larger than the system 
size, and the concept of local thermalization breaks down. 
Thus, the hydrodynamic evolution has to be stopped by 
applying the freeze-out criteria. We assume that freeze 
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out occurs at a fixed temperature Tp. We have used 
three different values of Tp (130, 140, and 150 MeV) 
in our calculation. By using the standard Cooper-Frye 
formalism [19], we calculate the invariant distribution of 
particles at the freeze-out hypersurface. In the Cooper- 
Frye formalism, the invariant distribution is given by the 
following equation: 



dN, 



9i 



(Pp dypTdprdcf) (27r)^ 



fi{p.u{x),x)p''d^af_, 



(12) 

where cJ'^cr^ is the outward normal vector on the freeze- 
out hypersurface S such that p^ fid^u^ is the local flux 
of particles with the momentum p through this surface, 
and the distribution function is 



1 



exp 



(E-tH(x))_ 



]±1 



(13) 



By applying Lorentz boost, we will get the value of lo- 
cal flow velocity u^{x) to the global reference frame by 
the substitution E— > p^u^^ where /ii(x) and T{x) are the 
chemical potential of particle species i and the local tem- 
perature along hypersurface S, respectively. From the 
invariant distribution, particle multiplicity (dNch/dy), 
mean pp (< pr >), elliptic flow (t'2), etc., can easily 
be computed. 



III. RESULTS 



Estimation of initial energy density in pp 
collisions at LHC 




FIG. 2: Variation of dNch/dy with impact parameter b, for 
^ = 14 TeV for C = 0.25 fm (dotted line) and ^ = 0.105 fm 
(solid line). 
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FIG. 3: (Color online) The minimum bias px spectra for three 
initial times {n = 0.2, 0.4, 0.6 fm), with I = 0.25 fm (dotted 
lines) and ^ = 0.105 fm (solid lines) at = 14 TeV. 
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FIG. 1: The filled circles are the RHIC data for the transverse 
momentum distribution for vr" in pp collisions at ^/s — 200 
GeV. The solid line is the hydrodynamical model fit to the 
data at diffuseness parameter ^ = 0.25 fm. 

In the absence of experimental guidance in pp colli- 
sions at the LHC {y/s = lA TeV), it is difficult to fix the 
initial energy density for hydrodynamical calculations 
in pp collisions at the LHC. We choose the initial energy 
density such that the experimentally measured or extrap- 
olated values of dNch/drj in pp collisions at given y^s is 



reproduced. We have fitted the experimentally measured 
Pt spectra of tt^ in pp collisions at y/s = 200 GeV from 
the STAR experiment at the RHIC, and they are shown 
in Figll] The filled circles are the RHIC data points for 
the transverse momentum distribution of tt~ . The solid 
line is the hydrodynamical model fit to the data with 
proton diffuseness parameter ^ = 0.25 fm. Note that the 
pp collisions are measured at minimum bias, without any 
centrality or impact-parameter dependence. To compare 
the experimental data with hydrodynamic simulations, 
we compute the minimum bias spectra as 



dN 

dyd^pj 



dN{b) 
dy<PpT 



db 



(14) 



where 



dN{b) 



is the TT invariant distribution at an im- 
We have integrated upto bmax = 1.6 
90% of thr INEL cross section. Gen- 



dyd^PT 

pact parameter b. 
fm, which covers ^ 
erally, one does not apply hydrodynamics in pp collisions. 
The system size possibly is too small for macroscopic con- 
cepts to be valid. However, in Ref. [l^l, there are some 
discussions on applications of hydrodynamical calcula- 
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be fixed only after detailed experimental investigations. 
Therefore, we perform our calculations for three different 
values of freeze-out temperature {Tp ~ 130, 140, and 150 
MeV) and three different values of initial time (tj = 0.2, 
0.4, and 0.6 fm) with two values of diffuseness parameter 
{£_ = 0.25 fm and ^ = 0.105 fm). 

The initial energy density (ej, the minimum bias av- 
erage multiplicity (^^^^) at y = 0, the meanpT, and the 
Pt integrated V2 in pp collisions at LHC energy (-y/s — 
14 TeV), are calculated and are summarized in Tables U 
and |lT] for various initial conditions at two values of dif- 
fuseness parameters. 



FIG. 4: (Color online) The minimum bias pr spectra for three 
freeze-out temperatures (Tf = 130, 140, 150 MeV) with ^ = 
0.25 fm (dotted lines) and ^ = 0.105 fm (solid lines) at \/s = 
14 TeV. 

tions to such systems. We observe that the pr spectra 
as measured by the STAR Collaboration at the RHIC 
{y/s ~ 200 GeV) for pp collisions is reasonably well ex- 
plained in the hydrodynamic model. In the low px re- 
gion < 0.5 GeV/c, the simulation underestimates the 
data because a large fraction of low px pions are prob- 
ably coming from various resonance decays, which are 
not considered in the simulation. Agreement with data 
in the low pT region will improve if resonance contribu- 
tions are accounted for. Hydrodynamical simulation also 
underestimates the data at highpr > 2 GeV/c, which in- 
dicates departure from ideal hydrodynamics because of 
dissipative effects. Also, at large pr, sources other than 
hydrodynamics contribute to particle production. 

Charged-particle pseudorapidity distribution in the 
central rapidity region in p-l-p and p+p interactions is 
expected to have a power-law dependence on centcr-of- 
mass energy [20l|. Therefore, we use the relation, 

— — =ax{y^i^Hc), (15) 

LHC 

where a and b are the parameters whose values are ob- 
tained as a = 0.7166 and b — 0.2171 by fitting the y/s de- 
pendence of INEL charged-particle pseudorapidity den- 
sity in the central rapidity region in p-l-p and p+p in- 
teractions [20j . The extrapolated minimum bias charged 
particles pseudorapidity density at y = at the LHC (-v/s 
= 14 TeV) is determined to be 5.69. For performing sim- 
ulations at the LHC, the initial energy density is fixed 
accordingly. The energy density required to reproduce 
the extrapolated dNch/drj is different for different values 
of initial time (t^) and freeze-out temperature {Tp). 

B. Mean multiplicity (dN/dy), mean pf {< pt >) 
and pt spectra 

The parameters that are used in this calculation for 
pp collisions need to be varied as these parameters can 



TABLE I: The minimum bias multiplicity dNch /dy, the mean 
Pt and the pr integrated «2 for various freeze-out tempera- 
tures for both diffuseness parameters at a fixed Ti. 



ri(fm) 


€(fm) 




(GcV//m3) 




<PT >(GeV/c) 


V2{%) 


0.6 


0.25 


130 


20.7 


5.68 


0.56 


0.34 


0.6 


0.25 


140 


26.2 


5.68 


0.57 


0.37 


0.6 


0.25 


150 


35.0 


5.68 


0.59 


0.40 


0.6 


0.105 


130 


31.4 


5.69 


0.68 


1.38 


0.6 


0.105 


140 


40.0 


5.69 


0.70 


1.38 


0.6 


0.105 


150 


53.8 


5.69 


0.73 


1.40 



TABLE II: The minimum bias multiplicity dNch/dy, the mean 
Pt and the pt integrated V2 for various initial times for both 
diffuseness parameters at a fixed Tf 



TF(MeV) 


« (fm) 


n (fm) 


(GeV//m3) 


dN^h/dy 


<PT >(GeV/c) 


V2{%) 


140 


0.25 


0.2 


99.7 


5.69 


0.60 


0.34 


140 


0.25 


0.4 


42.5 


5.68 


0.59 


0.35 


140 


0.25 


0.6 


26.2 


5.68 


0.57 


0.37 


140 


0.105 


0.2 


160.7 


5.69 


0.76 


1.32 


140 


0.105 


0.4 


66.5 


5.69 


0.72 


1.30 


140 


0.105 


0.6 


40.0 


5.69 


0.70 


1.38 



At a fixed initial time (r^ — 0.6 fm), when the freeze- 
out temperature is increased from 130 to 150 MeV, the 
mean pr is found to increase by 5.3% and 7.3% whereas 
the Pt integrated V2 is found to increase by 17.6% and 
1.4% for ^ = 0.25 and 0.105 fm, respectively. At a fixed 
freeze-out temperature (Tp = 140 MeV), when the initial 
time is increased from 0.2 to 0.6 fm, the mean px is found 
to decrease by 5.0% and 7.8%, whereas the pT integrated 
V2 is found to increase by 8.8% and 4.5% for ^ — 0.25 and 
0.105 fm, respectively. However, the meanpT is however, 
found to be consistently higher for the case with ^ = 0.105 
fm as compared to the case with ^ = 0.25 fm. The reason 
is understood. For the sharp proton density distribution 
(^ = 0.105 fm), fluid is initialized with higher energy 
density as compared to that when the surface is more 
diffused (^ — 0.25 fm). Increased transverse pressure in 
the fluid then leads to increased mean px- 

In Figl^J we have shown the impact-parameter depen- 
dence of dNch /dy at y = for both values of diffuseness 
parameter ^ = 0.105 fm and ^ = 0.25 fm. As expected. 
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FIG. 5: (Color online)The pr spectra for three different cen- 
tralities with ^ = 0.25 fm. 
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FIG. 6: (Color online)The pr spectra for three different cen- 
tralities with ^ — 0.105 fm. 



charged-particle rapidity density decreases with increas- 
ing impact parameters for both ^. This indicates that 
higher multipUcity events in pp coUisions are obtained at 
low-impact parameters, and a centrahty trigger based on 
multiphcity might select these interesting events. The 
rapidity density of charged particles at mid rapidity for 
b = is 13.79 for ^ = 0.105 fm and 12.01 for ^ = 0.25 fm. 
The impact-parameter dependence of rapidity density of 
a charged particle for two diffuseness parameters shows 
an interesting feature, that for ^ = 0.105 fm, the ra- 
pidity density of charged particles reduces faster as their 
configuration is close to the hard sphere; and above some 
impact parameter, the overlap region reduces faster. The 
result is according to our expectation. In peripheral col- 
lisions, the overlap region is comparatively large when 
the surface is more diffused. The impact parameter de- 
pendance of rapidity density of charged particles can be 
predicted to provide size information on the diffuseness 
parameter of the colliding protons. 

The minimum bias pr spectra of charged particles for 
different values of initial times and freeze-out tempera- 
ture are shown in FigsjS] ancQl respectively. Minimum 
bias pt spectra are obtained by integrating up to bmax 
= 1.6 fm, which, as previously indicated, covers ~90% 



of INEL cross section. It is observed that the px spectra 
are not affected much because of the change in freeze- 
out temperature and in initial time for a particular dif- 
fuseness parameter. However, the spectra are different 
for different diffuseness parameters. Note that, when- 
ever there is a change in either Tp at fixed or in at 
fixed Tp, the initial energy density (e^) is adjusted ac- 
cordingly so that the average multiplicity {dNch/dy) is 
always fixed. 




FIG. 7: Initial spatial eccentricity, Sx as a function of impact 
parameter for pp collisions at -^s — 14 TeV for ^ = 0.25 fm 
(dotted line) and ^ = 0.105 fm (solid line). 

As discussed earlier, the pseudorapidity density of the 
cgarged particle is 13.79 and 12.01 for ^ = 0.105 fm and ^ 
= 0.25 fm, respectively, at b = 0. This suggests that one 
can possibly determine two to three centralities experi- 
mentally on the basis of multiplicity. Therefore, we have 
studied the centrality (or impact parameter) dependence 
of spectra, V2 etc. for three centralities viz. to 20%, 
20 to 50% and 50 to 90%. The centrality dependence of 
Pt spectra for three centralities is shown in FigslS] ant^ 
for ^ = 0.25 fm and ^ — 0.105 fm, respectively. As in- 
dicated previously, we have neglected the resonance con- 
tribution. Resonance mainly contributes at the low px 
region. Centrality dependence of pt spectra is qualita- 
tively similar to that obtained in Au-|-Au collisions at 
the RHIC at both ^. As the collisions become more and 
more peripheral, the slope of the spectra gets steeper, 
which indicates reduced source temperature. The result 
is consistent with our expectation. Although the spectra 
looks similar for both ^, they differ in details (e.g., slope, 
total yield, etc.). 



C. Elliptic flow {v2) 

The particle azimuthal distribution can be constructed 
from different quantities such as transverse momentum, 
multiplicity, or transverse energy in relatively narrow 
(pseudo)-rapidity windows [2l| . For the non zero impact 
parameter, the azimuthal distribution of the particles in 
the reaction plane is not symmetric (anisotropic) in (p. 
We can decompose the distribution in terms of Fourier 
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expansion as 



E- 



d^p 2Trptdptdy 



{^(l + 2z;„cos(n0')}, (16) 



where (j)' is the angle of azimuth of the outgoing parti- 
cle with respect to the reaction plane. The coefficient 
of the second term in the Fourier expansion is known as 
the elliptic flow V2- The nonzero V2 describes the eccen- 
tricity of an ellipselike distribution. The origin of the 
measured anisotropy could be different: hydrodynami- 
cal flow, shadowing effect, both, etc. What they have in 
common is some collective behavior in the evolution of 
the multiparticle production process ^ . 

The initial spatial eccentricity of the reaction zone in 
the transverse plane changes with a change in the impact 
parameter. The V2 is sensitive to this change in spatial 
eccentricity and, consequently, to the change in the im- 
pact parameter. Initial spatial eccentricity is defined 
as 



(y2 + a;2)' 



(17) 



where the angular brackets denote energy density 
weighted averages at the initial time r^. In FiglTj ini- 
tial eccentricity Ex, for the two values of diffuseness pa- 
rameter ^ = 0.25 fm and ^ = 0.105 fm, are shown as 
a function of the impact parameter. For diffused pro- 
tons = 0.25 fm), the spatial eccentricity does not grow 
significantly. However, for the small diffuseness param- 
eter, ^ = 0.105 fm, the spatial eccentricity increases 
with the impact parameter and reaches a maximum at 
midcentral collisions. The behavior is similar to that of 
Au-I-Au collisions. 

The effect of varying initial conditions (t^ and Tp) on 
minimum bias differential elliptic flow (w2(pt)) is studied, 
and the results are depicted in Figs.|S]ancjni In Fig.|Sl we 
have plotted the minimum bias differential elliptic flow 
V2{pt) for three values of initial time (r^ = 0.2, 0.4, 0.6 
fm) for both ^. In Fig. ^ we have plotted the minimum 
bias differential elliptic flow V2{pt) for three values of 
freeze-out temperature {Tp = 130, 140, 150 MeV) for 
both ^. 

The minimum bias differential elliptic flow [v2{pt)], is 
found to be positive for both ^. The change in V2{pt) 
caused by the change in initial conditions (x; or Tp) is 
negligible. However, V2{pt) is found to be very large 
(reaches up to 24% &i px — i GeV/c) for the case with 
^ = 0.105 fm, whereas for ^ = 0.25 fm, it reaches a max- 
imum value up to 6% (at pt — i GeV/c). 

We have also studied the centrality dependence of pr 
integrated V2 and differential elliptic fiow V2{pt) for both 
^, for Ti = 0.6 fm and Tp = 140 MeV. The V2{pt) for 
three different centralities is shown in Fig. [lUl The dif- 
ferential elliptic flow V2{pt) is found to increase when 
we go toward mid peripheral, which reaches a maximum 
upto 33% (at pp = 3 GeV/c) for = 0.105 fm whereas 
it reaches a maximum value of 7% for ^ = 0.25 fm. The 



0.25 
0.2 

0.15 
0.1 

0.05 



- t,=0.105,Ti=0.6fm 




^ i,=0.105,T.i=0.4fm 




: l,=0.105,Xi=0.2fm . 




5=0.25, T,=0.6 fm 




i-0.25,i:-0.4 fm 




" i=0.25,Ti=0.2fm 




. , , 1 , . , . 1 , . , , 1 , , , . 1 . . , . 1 , , , 



0.5 



(4v) 



2.5 



FIG. 8: (Color online) The minimum bias differential elliptic 
flow V2{pt) for three initial times (r^ = 0.2, 0.4, 0.6 fm) for ^ 
= 0.25 fm (dotted lines) and ^ — 0.105 fm (solid lines). 



0.25 



0.2 



0.15 



0.1- 



0.05 



%=0.105, T ^=130 MeV 

- 5=0. 105, T^=140 MeV 

{,=0.105, T ^=150 MeV 
i,=0.25, 7^=130 MeV 
{,=0.25, T ^=140 MeV 
0.25, T =150 MeV 




1.5 2 

(GeV) 



FIG. 9: (Color online) The minimum bias differential elliptic 
flow V2{pt) for three freeze-out temperatures {Tp = 130, 140, 
150 MeV) for ^ = 0.25 fm (dotted lines) and ^ = 0.105 fm 
(solid lines). 



centrality dependence of W2(pt) is similar for both ^ apart 
from the differences in their respective values. 

The dependence of mean pp, integrated V2, and 
dNch/<fy with changing the centrality is summarized in 
Tables HH and |TV] for ^ = 0.25 fm and ^ = 0.105 fm, 
respectively. 



TABLE III: The minimum bias multiplicity dNch/dy, the 
mean pp, and the pr integrated V2 for three different cen- 
tralities at flxed Ti (0.6 fm), Tp (140 MeV) and ^ (0.25 fm). 



centrality 




<pp >(GeV/c) 


V2{Vo) 


0-20% 


10.78 


0.64 


0.13 


20-50% 


7.55 


0.60 


0.34 


50-90% 


3.45 


0.53 


0.29 



We observe that the pp integrated V2 shows centrality 
dependence as expected in a hydrodynamic model (i.e., 
increases when we go from central to midperipheral) , 
however, the value of pp integrated V2 is very small. For 
midcentral collisions, it reaches up to 1.34% and 0.34% 
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TABLE IV: The minimum bias multiplicity dNch/dy, the 
mean pr, and the pr integrated V2 for three different cen- 
tralities at fixed n (0.6 fm), Tf (140 MeV) and ^ (.105 fm). 



centrality 




<PT >(GeV/c) 


V2{%) 


0-20% 


11.74 


0.74 


0.66 


20-50% 


7.47 


0.71 


1.34 


50-90% 


3.26 


0.66 


1.23 
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FIG. 10: (Color online) The differential elliptic flow V2{pt) 
for three different centralities for ^ — 0.25 fm (dotted line) 
and ^ — 0.105 fm (solid line). 

for ^ — 0.105 fm and ^ — 0.25 fm, respectively. 

IV. SUMMARY AND CONCLUSION 

To summarize, we have presented a hydrodynamical 
model study for pp collisions at LHC {y/s — 14 TeV) en- 



ergy. The initial energy density for the hydrodynamic 
calculation is obtained by requiring a condition that the 
extrapolated rapidity density of charged particles should 
be reproduced in the model. After studying the basic 
properties of particle production (e.g., spectra, dN/dy, 
(Pt), and their centrality dependence), we have extracted 
V2 and its dependence on pt and centralities. As the ini- 
tial spatial eccentricity depends strongly on the shape of 
partonic density inside the proton, we have studied two 
cases of surface diffuseness parameters in Woods-Saxon 
profile. We have studied the effect of changing initial time 
Ti and freeze-out temperature Tp on elliptic flow. It is 
found that the effect is negligible. The minimum bias 
differential elliptic flow V2{pt) is found to be very large 
and reaches up to 24% for the case with diffuseness pa- 
rameter, ^ = 0.105 fm. When the diffuseness parameter 
is 0.25 fm, the minimum bias differential elliptic flow is 
found to reach up to 6%. At LHC energy, the energy 
density and particle multiplicity are expected to be sub- 
stantially large, which suggests the creation of a system 
that shows collective properties. First results from the 
LHC 2^ have already shown charged multiplicity at 
the midrapidity region going up to 60. Therefore, we 
suggest to measure V2 and its centrality dependence in 
pp collisions at the LHC. The centrality dependence can 
be studied by selecting events with different multiplicity 
bins. The estimation of V2 at LHC energy might be af- 
fected by the presence of nonflow effects (e.g. jets), but 
methods based on a scalar product, a cumulant or a stan- 
dard event plane with a large jy-gap can be applied. The 
centrality dependence of V2 can be used to obtain the 
surface diffuseness parameter of protons in pp collisions. 
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